Non-canonical functions of UHRF1 maintain DNA methylation homeostasis in cancer cells.

DNA methylation is an essential epigenetic chromatin modification, and its maintenance in mammals requires the protein UHRF1. It is yet unclear if UHRF1 functions solely by stimulating DNA methylation maintenance by DNMT1, or if it has important additional functions. Using degron alleles, we show that UHRF1 depletion causes a much greater loss of DNA methylation than DNMT1 depletion. This is not caused by passive demethylation as UHRF1-depleted cells proliferate more slowly than DNMT1-depleted cells. Instead, bioinformatics, proteomics and genetics experiments establish that UHRF1, besides activating DNMT1, interacts with DNMT3A and DNMT3B and promotes their activity. In addition, we show that UHRF1 antagonizes active DNA demethylation by TET2. Therefore, UHRF1 has non-canonical roles that contribute importantly to DNA methylation homeostasis; these findings have practical implications for epigenetics in health and disease.

Transcription factor-mediated direct cellular reprogramming yields cell-type specific DNA methylation signature.

Direct reprogramming, inducing the conversion of one type of somatic cell into another by the forced expression of defined transcription factors, is a technology with anticipated medical applications. However, due to the many unresolved aspects of the induction mechanisms, it is essential to thoroughly analyze the epigenomic state of the generated cells. Here, we performed comparative genome-wide DNA methylation analyses of mouse embryonic fibroblasts (MEFs) and cells composing organoids formed by intestinal stem cells (ISCs) or induced ISCs (iISCs) that were directly induced from MEFs. We found that the CpG methylation state was similar between cells forming ISC organoids and iISC organoids, while they differed widely from those in MEFs. Moreover, genomic regions that were differentially methylated between ISC organoid- and iISC organoid-forming cells did not significantly affect gene expression. These results demonstrate the accuracy and safety of iISC induction, leading to the medical applications of this technology.

DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins in liver and skeletal muscle.

Each tissue has a dominant set of functional proteins required to mediate tissue-specific functions. Epigenetic modifications, transcription, and translational efficiency control tissue-dominant protein production. However, the coordination of these regulatory mechanisms to achieve such tissue-specific protein production remains unclear. Here, we analyzed the DNA methylome, transcriptome, and proteome in mouse liver and skeletal muscle. We found that DNA hypomethylation at promoter regions is globally associated with liver-dominant or skeletal muscle-dominant functional protein production within each tissue, as well as with genes encoding proteins involved in ubiquitous functions in both tissues. Thus, genes encoding liver-dominant proteins, such as those involved in glycolysis or gluconeogenesis, the urea cycle, complement and coagulation systems, enzymes of tryptophan metabolism, and cytochrome P450-related metabolism, were hypomethylated in the liver, whereas those encoding-skeletal muscle-dominant proteins, such as those involved in sarcomere organization, were hypomethylated in the skeletal muscle. Thus, DNA hypomethylation characterizes genes encoding tissue-dominant functional proteins.

The Epigenetic Legacy of Maternal Protein Restriction: Renal Ptger1 DNA Methylation Changes in Hypertensive Rat Offspring.

Nutrient imbalances during gestation are a risk factor for hypertension in offspring. Although the effects of prenatal nutritional deficiency on the development of hypertension and cardiovascular diseases in adulthood have been extensively documented, its underlying mechanisms remain poorly understood. In this study, we aimed to elucidate the precise role and functional significance of epigenetic modifications in the pathogenesis of hypertension. To this end, we integrated methylome and transcriptome data to identify potential salt-sensitive hypertension genes using the kidneys of stroke-prone spontaneously hypertensive rat (SHRSP) pups exposed to a low-protein diet throughout their fetal life. Maternal protein restriction during gestation led to a positive correlation between DNA hypermethylation of the renal prostaglandin E receptor 1 (Ptger1) CpG island and high mRNA expression of Ptger1 in offspring, which is consistently conserved. Furthermore, post-weaning low-protein or high-protein diets modified the Ptger1 DNA hypermethylation caused by fetal malnutrition. Here, we show that this epigenetic variation in Ptger1 is linked to disease susceptibility established during fetal stages and could be reprogrammed by manipulating the postnatal diet. Thus, our findings clarify the developmental origins connecting the maternal nutritional environment and potential epigenetic biomarkers for offspring hypertension. These findings shed light on hypertension prevention and prospective therapeutic strategies.

A highly efficient scheme for library preparation from single-stranded DNA.

Although methods for sequencing library preparation from double-stranded DNA are well established, those from single-stranded DNA (ssDNA) have not been well studied. Further, the existing methods have limitations in efficiency and yield. Therefore, we developed a highly efficient procedure for sequencing library preparation from ssDNA. In this method, the first adaptor tagging of ssDNA is performed using terminal deoxyribonucleotidyl transferase (TdT)-assisted adenylate connector-mediated ssDNA (TACS) ligation, which we reported recently. After complementary strand synthesis using the adaptor-tagged ssDNA, second adaptor tagging via Vaccinia virus topoisomerase I (VTopoI or TOPO)-based adaptor ligation is performed. With additional steps for degradation, repression, and removal of the adaptor dimer, the proposed TACS-TOPO scheme realizes adaptor dimer-free sequencing library preparation from ssDNA samples of 24 pg. The TACS-TOPO scheme was successfully applied to cell-free DNA analysis with amplification-free library preparation from 50 µL of human serum. A modified TACS-TOPO scheme was also applied to DNA extracted from ancient human bones, bringing two to eight times more library yields than those using a conventional library preparation protocol. The procedures for preparing VTopoI and its complex with a double-stranded oligonucleotide adaptor are also described. Overall, the proposed TACS-TOPO scheme can facilitate practical and sensitive sequencing analysis of ssDNA.

A genetic screen identifies BEND3 as a regulator of bivalent gene expression and global DNA methylation.

Epigenetic mechanisms are essential to establish and safeguard cellular identities in mammals. They dynamically regulate the expression of genes, transposable elements and higher-order chromatin structures. Consequently, these chromatin marks are indispensable for mammalian development and alterations often lead to disease, such as cancer. Bivalent promoters are especially important during differentiation and development. Here we used a genetic screen to identify new regulators of a bivalent repressed gene. We identify BEND3 as a regulator of hundreds of bivalent promoters, some of which it represses, and some of which it activates. We show that BEND3 is recruited to a CpG-containg consensus site that is present in multiple copies in many bivalent promoters. Besides having direct effect on the promoters it binds, the loss of BEND3 leads to genome-wide gains of DNA methylation, which are especially marked at regions normally protected by the TET enzymes. DNA hydroxymethylation is reduced in Bend3 mutant cells, possibly as consequence of altered gene expression leading to diminished alpha-ketoglutarate production, thus lowering TET activity. Our results clarify the direct and indirect roles of an important chromatin regulator, BEND3, and, more broadly, they shed light on the regulation of bivalent promoters.

A genome-wide screen reveals new regulators of the 2-cell-like cell state.

In mammals, only the zygote and blastomeres of the early embryo are totipotent. This totipotency is mirrored in vitro by mouse '2-cell-like cells' (2CLCs), which appear at low frequency in cultures of embryonic stem cells (ESCs). Because totipotency is not completely understood, we carried out a genome-wide CRISPR knockout screen in mouse ESCs, searching for mutants that reactivate the expression of Dazl, a gene expressed in 2CLCs. Here we report the identification of four mutants that reactivate Dazl and a broader 2-cell-like signature: the E3 ubiquitin ligase adaptor SPOP, the Zinc-Finger transcription factor ZBTB14, MCM3AP, a component of the RNA processing complex TREX-2, and the lysine demethylase KDM5C. All four factors function upstream of DPPA2 and DUX, but not via p53. In addition, SPOP binds DPPA2, and KDM5C interacts with ncPRC1.6 and inhibits 2CLC gene expression in a catalytic-independent manner. These results extend our knowledge of totipotency, a key phase of organismal life.

Crucial role of iron in epigenetic rewriting during adipocyte differentiation mediated by JMJD1A and TET2 activity.

Iron metabolism is closely associated with the pathogenesis of obesity. However, the mechanism of the iron-dependent regulation of adipocyte differentiation remains unclear. Here, we show that iron is essential for rewriting of epigenetic marks during adipocyte differentiation. Iron supply through lysosome-mediated ferritinophagy was found to be crucial during the early stage of adipocyte differentiation, and iron deficiency during this period suppressed subsequent terminal differentiation. This was associated with demethylation of both repressive histone marks and DNA in the genomic regions of adipocyte differentiation-associated genes, 　including Pparg, which encodes PPARγ, the master regulator of adipocyte differentiation. In addition, we identified several epigenetic demethylases to be responsible for iron-dependent adipocyte differentiation, with the histone demethylase jumonji domain-containing 1A and the DNA demethylase ten-eleven translocation 2 as the major enzymes. The interrelationship between repressive histone marks and DNA methylation was indicated by an integrated genome-wide association analysis, and was also supported by the findings that both histone and DNA demethylation were suppressed by either the inhibition of lysosomal ferritin flux or the knockdown of iron chaperone poly(rC)-binding protein 2. In summary, epigenetic regulations through iron-dependent control of epigenetic enzyme activities play an important role in the organized gene expression mechanisms of adipogenesis.

Conservation and divergence of canonical and non-canonical imprinting in murids.

BACKGROUND: Genomic imprinting affects gene expression in a parent-of-origin manner and has a profound impact on complex traits including growth and behavior. While the rat is widely used to model human pathophysiology, few imprinted genes have been identified in this murid. To systematically identify imprinted genes and genomic imprints in the rat, we use low input methods for genome-wide analyses of gene expression and DNA methylation to profile embryonic and extraembryonic tissues at allele-specific resolution. RESULTS: We identify 14 and 26 imprinted genes in these tissues, respectively, with 10 of these genes imprinted in both tissues. Comparative analyses with mouse reveal that orthologous imprinted gene expression and associated canonical DNA methylation imprints are conserved in the embryo proper of the Muridae family. However, only 3 paternally expressed imprinted genes are conserved in the extraembryonic tissue of murids, all of which are associated with non-canonical H3K27me3 imprints. The discovery of 8 novel non-canonical imprinted genes unique to the rat is consistent with more rapid evolution of extraembryonic imprinting. Meta-analysis of novel imprinted genes reveals multiple mechanisms by which species-specific imprinted expression may be established, including H3K27me3 deposition in the oocyte, the appearance of ZFP57 binding motifs, and the insertion of endogenous retroviral promoters. CONCLUSIONS: In summary, we provide an expanded list of imprinted loci in the rat, reveal the extent of conservation of imprinted gene expression, and identify potential mechanisms responsible for the evolution of species-specific imprinting.

Post-bisulfite Adaptor Tagging with a Highly Efficient Single-Stranded DNA Ligation Technique.

Post-bisulfite adaptor tagging (PBAT) is a procedure for efficiently preparing a sequencing library for whole-genome bisulfite sequencing (WGBS). The original version of the PBAT protocol was highly efficient, such that it helped realize library preparation from samples of limited amounts. However, two rounds of random priming reactions employed in the original protocol limited further improvement of the PBAT protocol in terms of read length and mapping rate. In this chapter, an improved version of the PBAT protocol called tPBAT is described.

Post-bisulfite Adaptor Tagging Based on an ssDNA Ligation Technique (tPBAT).

Post-bisulfite adaptor tagging (PBAT) is a concept that enables the preparation of an efficient sequencing library from bisulfite-treated DNA, and it also means the protocol implemented the concept. Although the previous PBAT or rPBAT was sensitive enough for single-cell methylome analysis, the protocol had several drawbacks owing to the repeated random priming reactions. To resolve these problems, we developed a unique single-strand DNA ligation technique, termed TACS ligation, and established a new protocol called tPBAT. With tPBAT, the data quality improved, with a longer insert and higher mapping rate than that obtained with rPBAT. In addition, paired-end sequencing and indexing were supported by the default. In this chapter, the tPBAT protocol is introduced, and a thorough description of its application to small samples is provided.

Persistent epigenetic alterations in transcription factors after a sustained virological response in hepatocellular carcinoma.

Background and Aim: The risk of hepatocellular carcinoma (HCC) persists in a condition of sustained virologic response (SVR) after hepatitis C virus (HCV) eradication. Comprehensive molecular analyses were performed to test the hypothesis that epigenetic abnormalities present after an SVR play a role in hepatocarcinogenesis. Methods: Whole-genome methylome and RNA sequencing were performed on HCV, SVR, and healthy liver tissue. Integrated analysis of the sequencing data focused on expression changes in transcription factors and their target genes, commonly found in HCV and SVR. Identified expression changes were validated in demethylated cultured HCC cell lines and an independent validation cohort. Results: The coincidence rates of the differentially methylated regions between the HCV and SVR groups were 91% in the hypomethylated and 71% in the hypermethylated regions in tumorous tissues, and 37% in the hypomethylated and 36% in the hypermethylated regions in non-tumorous tissues. These results indicate that many epigenomic abnormalities persist even after an SVR was achieved. Integrated analysis identified 61 transcription factors and 379 other genes that had methylation abnormalities and gene expression changes in both groups. Validation cohort specified gene expression changes for 14 genes, and gene ontology pathway analysis revealed apoptotic signaling and inflammatory response were associated with these genes. Conclusion: This study demonstrates that DNA methylation abnormalities, retained after HCV eradication, affect the expression of transcription factors and their target genes. These findings suggest that DNA methylation in SVR patients may be functionally important in carcinogenesis, and could serve as biomarkers to predict HCC occurrence.

Identification, expression, and purification of DNA cytosine 5-methyltransferases with short recognition sequences.

BACKGROUND: DNA methyltransferases (MTases) are enzymes that induce methylation, one of the representative epigenetic modifications of DNA, and are also useful tools for analyzing epigenomes. However, regarding DNA cytosine 5-methylation, MTases identified so far have drawbacks in that their recognition sequences overlap with those for intrinsic DNA methylation in mammalian cells and/or that the recognition sequence is too long for fine epigenetic mapping. To identify MTases with short recognition sequences that never overlap with the CG dinucleotide, we systematically investigated the 25 candidate enzymes identified using a database search, which showed high similarity to known cytosine 5-MTases recognizing short sequences. RESULTS: We identified MTases with six new recognition sequences, including TCTG, CC, CNG, TCG, GCY, and GGCA. Because the recognition sequence never overlapped with the CG dinucleotide, MTases recognizing the CC dinucleotide were promising. CONCLUSIONS: In the current study, we established a procedure for producing active CC-methylating MTases and applied it to nucleosome occupancy and methylome sequencing to prove the usefulness of the enzyme for fine epigenetic mapping. MTases that never overlap with CG dinucleotides would allow us to profile multiple epigenomes simultaneously.

An efficient procedure for the recovery of DNA from formalin-fixed paraffin-embedded tissue sections.

With the advent of new molecular diagnostic techniques, retrieving DNA from the formalin-fixed paraffin-embedded (FFPE) tissues has become an essential yet challenging step for efficient downstream processes. Owing to low quality and quantity of DNA retrieved from the FFPE sections, the process is often impractical and needs significant improvements. Here, we established an efficient method for the purification of DNA from FFPE specimens by optimizing incubation temperature, incubation time, and the concentration of a formalin scavenger tris(hydroxymethyl)aminomethane (Tris) for reverse-crosslinking. The optimized method, named "Highly concentrated Tris-mediated DNA extraction" (HiTE), yielded three times the DNA yield per tissue slice compared with a representative DNA extraction kit. Moreover, the use of HiTE-extracted DNA increased the yield of the sequencing library three times and accordingly yielded a log higher and more reproducible sequencing library compared with that obtained using the commonly used commercial kit. The sequencing library prepared from HiTE-extracted FFPE-DNA had longer inserts and produced reads that evenly covered the reference genome. Successful application of HiTE-extracted FFPE-DNA for whole-genome and targeted gene panel sequencing indicates its practical usability.

Clonal germinal center B cells function as a niche for T-cell lymphoma.

Angioimmunoblastic T-cell lymphoma (AITL) is proposed to be initiated by age-related clonal hematopoiesis (ACH) with TET2 mutations, whereas the G17V RHOA mutation in immature cells with TET2 mutations promotes the development of T follicular helper (TFH)-like tumor cells. Here, we investigated the mechanism by which TET2-mutant immune cells enable AITL development using mouse models and human samples. Among the 2 mouse models, mice lacking Tet2 in all the blood cells (Mx-Cre × Tet2flox/flox × G17V RHOA transgenic mice) spontaneously developed AITL for approximately up to a year, while mice lacking Tet2 only in the T cells (Cd4-Cre × Tet2flox/flox × G17V RHOA transgenic mice) did not. Therefore, Tet2-deficient immune cells function as a niche for AITL development. Single-cell RNA-sequencing (scRNA-seq) of >50 000 cells from mouse and human AITL samples revealed significant expansion of aberrant B cells, exhibiting properties of activating light zone (LZ)-like and proliferative dark zone (DZ)-like germinal center B (GCB) cells. The GCB cells in AITL clonally evolved with recurrent mutations in genes related to core histones. In silico network analysis using scRNA-seq data identified Cd40-Cd40lg as a possible mediator of GCB and tumor cell cluster interactions. Treatment of AITL model mice with anti-Cd40lg inhibitory antibody prolonged survival. The genes expressed in aberrantly expanded GCB cells in murine tumors were also broadly expressed in the B-lineage cells of TET2-mutant human AITL. Therefore, ACH-derived GCB cells could undergo independent clonal evolution and support the tumorigenesis in AITL via the CD40-CD40LG axis.

A GCDGC-specific DNA (cytosine-5) methyltransferase that methylates the GCWGC sequence on both strands and the GCSGC sequence on one strand.

5-Methylcytosine is one of the major epigenetic marks of DNA in living organisms. Some bacterial species possess DNA methyltransferases that modify cytosines on both strands to produce fully-methylated sites or on either strand to produce hemi-methylated sites. In this study, we characterized a DNA methyltransferase that produces two sequences with different methylation patterns: one methylated on both strands and another on one strand. M.BatI is the orphan DNA methyltransferase of Bacillus anthracis coded in one of the prophages on the chromosome. Analysis of M.BatI modified DNA by bisulfite sequencing revealed that the enzyme methylates the first cytosine in sequences of 5'-GCAGC-3', 5'-GCTGC-3', and 5'-GCGGC-3', but not of 5'-GCCGC-3'. This resulted in the production of fully-methylated 5'-GCWGC-3' and hemi-methylated 5'-GCSGC-3'. M.BatI also showed toxicity when expressed in E. coli, which was caused by a mechanism other than DNA modification activity. Homologs of M.BatI were found in other Bacillus species on different prophage like regions, suggesting the spread of the gene by several different phages. The discovery of the DNA methyltransferase with unique modification target specificity suggested unrevealed diversity of target sequences of bacterial cytosine DNA methyltransferase.

ChIP-Atlas 2021 update: a data-mining suite for exploring epigenomic landscapes by fully integrating ChIP-seq, ATAC-seq and Bisulfite-seq data.

ChIP-Atlas (https://chip-atlas.org) is a web service providing both GUI- and API-based data-mining tools to reveal the architecture of the transcription regulatory landscape. ChIP-Atlas is powered by comprehensively integrating all data sets from high-throughput ChIP-seq and DNase-seq, a method for profiling chromatin regions accessible to DNase. In this update, we further collected all the ATAC-seq and whole-genome bisulfite-seq data for six model organisms (human, mouse, rat, fruit fly, nematode, and budding yeast) with the latest genome assemblies. These together with ChIP-seq data can be visualized with the Peak Browser tool and a genome browser to explore the epigenomic landscape of a query genomic locus, such as its chromatin accessibility, DNA methylation status, and protein-genome interactions. This epigenomic landscape can also be characterized for multiple genes and genomic loci by querying with the Enrichment Analysis tool, which, for example, revealed that inflammatory bowel disease-associated SNPs are the most significantly hypo-methylated in neutrophils. Therefore, ChIP-Atlas provides a panoramic view of the whole epigenomic landscape. All datasets are free to download via either a simple button on the web page or an API.

Pan-cancer methylome analysis for cancer diagnosis and classification of cancer cell of origin.

The accurate and early diagnosis and classification of cancer origin from either tissue or liquid biopsy is crucial for selecting the appropriate treatment and reducing cancer-related mortality. Here, we established the CAncer Cell-of-Origin (CACO) methylation panel using the methylation data of the 28 types of cancer in The Cancer Genome Atlas (7950 patients and 707 normal controls) as well as healthy whole blood samples (95 subjects). We showed that the CACO methylation panel had high diagnostic potential with high sensitivity and specificity in the discovery (maximum AUC = 0.998) and validation (maximum AUC = 1.000) cohorts. Moreover, we confirmed that the CACO methylation panel could identify the cancer cell type of origin using the methylation profile from liquid as well as tissue biopsy, including primary, metastatic, and multiregional cancer samples and cancer of unknown primary, independent of the methylation analysis platform and specimen preparation method. Together, the CACO methylation panel can be a powerful tool for the classification and diagnosis of cancer.

Neuronal activation modulates enhancer activity of genes for excitatory synaptogenesis through de novo DNA methylation.

Post-mitotic neurons do exhibit DNA methylation changes, contrary to the longstanding belief that the epigenetic pattern in terminally differentiated cells is essentially unchanged. While the mechanism and physiological significance of DNA demethylation in neurons have been extensively elucidated, the occurrence of de novo DNA methylation and its impacts have been much less investigated. In the present study, we showed that neuronal activation induces de novo DNA methylation at enhancer regions, which can repress target genes in primary cultured hippocampal neurons. The functional significance of this de novo DNA methylation was underpinned by the demonstration that inhibition of DNA methyltransferase (DNMT) activity decreased neuronal activity-induced excitatory synaptogenesis. Overexpression of WW and C2 domain-containing 1 (Wwc1), a representative target gene of de novo DNA methylation, could phenocopy this DNMT inhibition-induced decrease in synaptogenesis. We found that both DNMT1 and DNMT3a were required for neuronal activity-induced de novo DNA methylation of the Wwc1 enhancer. Taken together, we concluded that neuronal activity-induced de novo DNA methylation that affects gene expression has an impact on neuronal physiology that is comparable to that of DNA demethylation. Since the different requirements of DNMTs for germ cell and embryonic development are known, our findings also have considerable implications for future studies on epigenomics in the field of reproductive biology.